Correlation-function asymptotic expansions: pre-factors universality 
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We show that the pre-factors of all terms of the one-dimensional Hubbard model correlation- 
function asymptotic expansions have an universal form, as the corresponding critical exponents. In 
addition to calculating such pre-factors, our study clarifies the relation of the low-energy Tomonaga- 
Luttinger liquid behavior to the scattering mechanisms which control the spectral properties of the 
model at all energy scales. Our results are of general nature for many integrable interacting models 
and provide a broader understanding of the unusual properties of quasi-lD nanostructures, organic 
conductors, and optical lattices of fermionic atoms. 
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Over the past 25 years it has been found that 
the low-energy physics of a variety of models of one- 
dimensional (ID) correlated electrons can be described 
by the Tomonaga-Luttinger liquid (TLL) theory . Im- 
portantly, the low-energy TLL universal behavior was ob- 
served in different real materials and systems, as for in- 
stance in carbon nanotubes 0, ballistic wires quasi- 
lD organic conductors 0, and quasi-lD quantum gases 
of ultracold fermionic atoms j^. On the other hand, the 
low-energy phases of some quasi-lD compounds are not 
metallic and correspond to broken-symmetry states [^. 
Recently, the resolution of photoemission experiments 
has improved, and the normal state of these compounds 
was found to display exotic spectral properties ^] . How- 
ever, such a metallic phase refers to finite energies and is 
not described by the TLL theory. 

The ID Hubbard model is one of the few realistic mod- 
els for correlated electrons in a discrete lattice for which 
one can exactly calculate all the energy eigenstates and 
their energies [g]- It includes a first-neighbor transfer- 
integral for electron hopping along the chain and an ef- 
fective on-site repulsion U. For finite energy values, the 
metallic phase of this model is beyond a TLL description 
and thus the study of spectral functions is a very involved 
many-electron problem. Fortunately, the recently intro- 
duced pseudofermion dynamical theory (PDT) provides 
explicit expressions for these functions Q. Moreover, 
the theory describes successfully the unusual spectral fea- 
tures of quasi-lD compounds for the whole finite-energy 
band width j^]. Recently, consistent results were ob- 
tained by numerical techniques 0|. Furthermore, when 
combined with the Renormalization Group, the use of 
the PDT reveals that a system of weakly coupled Hub- 
bard chains is suitable for the successful description of 
the phase diagram observed in quasi-lD doped Mott- 
Hubbard insulators 11]. In turn, the low-energy physics 
of the model corresponds to the universal TLL behavior 
and was studied by different techniques, such as bosoniza- 
tion and conformal-field theory (CFT) ^1^]. 



There are many investigations where the low-energy 
conformal invariance was combined with the model ex- 
act Bethe-ansatz solution in the study of the asymp- 
totics of correlation functions and related quantities 
Q El E 113, El El- However, the relation of the 
low-energy TLL behavior to the microscopic scattering 
mechanisms which control the unusual spectral proper- 
ties of the model at all energy scales remains an inter- 
esting open problem, which we address in this paper. 
While CFT and bosonization techniques do not provide 
correlation-function expressions for finite energy, here we 
use the general PDT to derive the asymptotic expansions 
for the model correlation functions. Such asymptotic ex- 
pansions go beyond those provided by other methods, 
once we can derive explicit expressions for the pre-factors 
of all the expansion terms. Moreover, we find that such 
pre-factors have an universal form, as the corresponding 
critical exponents. Our results describe the emergence 
of the TLL low-energy physics in terms of the general 
non-perturbative microscopic scattering mechanisms of 
the model at all energy scales and lead to a broader un- 
derstanding of the unusual properties observed in low- 
dimensional materials and nanostructures H Q and 
systems of interacting ultracold fermionic atoms in ID 
optical lattices 

In our study we consider the ID Hubbard model with 
periodic boundary conditions, number of lattice sites Na 
large and even, and units such that the transfer integral 
and the Planck and electronic-lattice constants are one. 
Thus, the lattice length is L = Na- Let us consider the 
general correlation function. 
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where Bj^{k,uj) is the A/'-electron spectral function, 

B'^{k,Lu)^Y. \{f\0'Ak)\GS)\'6ii.~lAEf). (2) 
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the general A/'-electron operators Oj^^{k) = 0\j-{k) and 
Ojf'ik) = Ojsf{k) carry momentum k, the / summation 
runs over the excited energy eigenstates, and A_Ey = 
[Ef—Ecs] where the energy Ef corresponds to these final 
states and Eqs is the ground-state energy. Most common 
examples are the operator Oi(fc) = Ck^a and different 
choices of charge, spin, and Cooper-pair A/" = 2 operators. 
We consider that the electronic density n = rif -f and 
spin density m = — ni are in the range < n < 1 and 
< m < n and thus the Fermi momenta are given by 
kp — Tm/2 and kp„ = nn^. 

The Fourier transform xVl'^' '^^ function Q rel- 
ative to both the momentum k and energy lo can be ex- 
pressed in terms of the corresponding Fourier transform 
Bj^{x,t) of the spectral function (|2Jl as follows, 

x'^ix,t) = -i2neilt)B'^{x,t), (3) 

where here and in other expressions of this paper 9{y) = 1 
for y > and 9{y) = for y < 0. 

Our main goal is the study of the general correlation 
function (O asymptotic expansion. Such an asymptotic 
expansion is controlled by the low-energy behavior of the 
function Q . From the studies of Ref. [22| one knows that 
the excitation spectrum is such that for the densities con- 
sidered here, the excited energy eigenstates which span 
the low-energy subspace can be represented by occupancy 
configurations of the c and si pseudoparticles, which are 
defined in terms of the electrons in that reference. We 
denote the si branch by s and use a general index a 
such that a — c, s. The a pseudoparticles have discrete 
bare- momentum values qj such that qj+i — qj — 2tt/L 
and energy residual-interaction terms associated with / 
functions fa,a'{q,q')- In the continuum limit the bare- 
momentum values q exist in the range \q\ < q^ where 
except for 1/L corrections gj? and read q^ — tt and 
5° = kp'i and play the role of Brillouin zone limits. It is 
useful to consider the index t = ±1 associated with the 
right (i = +1) and left = a, t Fermi points, which 
except for 1/L corrections are given by iq^p^ = i2kp and 
Lq^p^ = ikpi for the ground state. That state corresponds 
to a-band bare-momentum densely packed occupancy 
configurations such that |g| < qp^ and qp^ < \q\ < g° 
for a pseudofcrmions and a pseudofermion holes, respec- 
tively. (All a and t sums or products appearing in the 
expressions provided below run over the values a = c, s 
and t = +1, — 1, respectively.) 

The pseudoparticle - pseudofermion unitary transfor- 
mation plays a key role in the PDT. For the low-energy 
subspace, such a transformation maps the a pseudopar- 
ticles or pseudoparticle holes onto a pseudofcrmions or 
pseudofermion holes, by introducing shifts, Q%{qj)/L, 
of order 1/L in the discrete bare- momentum values 
and leaving all other pseudoparticle properties invariant 

H^l . The a pseudofermions have discrete canonical- 
momentum values qj = qj -\- Q'^{qj)/L and no energy 



residual-interaction terms. The point is that while for 
the electrons the present problem is strongly correlated, 
for the pseudofermions the shift Q^{q)/L results from 
zero-momentum forward scattering only. For the ground 
state, qj — qj and the pseudoparticle and pseudofermion 
are the same object. Each ground-state - excited-state 
transition involves a scattering event whose scattering 
centers are the a' , g' pseudofermions and pseudofermion 
holes created under the transition and the scatterers are 
all a,q pseudofermions and pseudofermion holes of the 
excited state, which acquire an overall scattering phase 
shift Q*(g)/2 = Eo,E^,7r$„„-((7,(z')AiV„,(g'). Here 
ANa'{q') is the excited-state a'-band bare- momentum 
distribution function deviation and ±7r (g, g') is 
a two-pseudofermion phase shift acquired by the a, q 
pseudofermion or hole scatterers for each a',q' pseud- 
ofermion or hole (— ) scattering center created under 
the transition. The value of such a phase shift is in the 
range — 7r/2 < tt (g, g') < +tt/2 and its U, n, and m 
dependence is uniquely defined by solution of a system 
of integral equations. The overall phase shift reads, 

Q„(g)/2 = Q0/2 + Q*(g)/2, (4) 

where Q^/2 — 0, ±7r/2 is a scattering-less contribution 
which has a single and well-defined value for the whole 
low-energy excitation subspace of each correlation func- 
tion. Such a subspace is spanned by all energy eigenstates 
with the same value for the number deviations ANc — 
AN and A7V,. = AN^ such that AN^' = J2q' ^^a' W) = 
X;,AAf^ ^. In contrast to tt $„„/(g, g'), Q*(g)/2 is a 
functional whose value depends on the excited state 
trough the deviations AiVa/(g'). For the present low- 
energy problem, such deviations refer to g' values in the 
vicinity of the Fermi points only and thus the number 
AiVf,^ = AiV";^ -f iQO,/27r is such that AN°f^ is the 
deviation in the number of a' pseudofermions at the a', l 
Fermi point. The a'-branch current-number deviation 
reads AJ^, = ^J2l '•^-^a' l- The low-energy subspace 
contains several J-suhspaces, which differ at least in one 
of the two values {Ajf,AJf}. The J-ground state is 
the lowest-energy eigenstate of a J-subspace. Its bare- 
momentum occupancy configurations are densely packed 
with Fermi points qpa,L = t-qpa + Aqpa^t- The deviation 
Aqpa^L and the excitation momentum of that state read, 

Ag^„, = .| + .Ajf ) ; fc^^ = ^ 2gL^Ajf, . 

a' 

(5) 

The momentum kg is generated by zero-energy and 
finite-momentum elementary processes, which create 
a' pseudofermions (AA^°;^ > 0) or a' pseudofermion 

holes (AAf°;^ < 0) at least at one of the four a',L 
Fermi points. The corresponding J-subspace is spanned 
by energy eigenstates generated from the J-ground 
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state by small- momentum and low-energy "particle- 
hole" pseudofermion processes in the vicinity of the 
a' , L Fermi points. Such processes conserve the set 
of {ANc, ANs, AJ[ , AJf} deviation values. A crucial 
point for the low-energy scattering properties is that 
the a' pseudofermions and holes created by such small- 
momentum and low-energy processes are not active scat- 
tering centers, once the phase shifts generated by the 
created pseudofermions exactly cancel those originated 
by creation of the corresponding holes. It follows that 
the overall phase shift has for each a pseudofermion 
or hole scatterer of bare-momentum q the the same value 
for all excited states spanning the same J-subspace with. 



(6) 

A mechanism which plays a central role in the spec- 
tral properties at all energy scales is the occurrence of 
a shift in the value of the four a, l canonical-momentum 
Fermi-points, as a result of each ground-state - excited- 
J-ground-state transition. The square of that shift reads, 



2A'^ = (AqFo.^J 



,2tt 



= [lAN. 



O.F 



27r 



(7) 



At low energy, consistently with the form of the scatter- 
ing phase shift ©, 2A^ simplifies to. 



2A^ = 2A^(AA^„A7V„Ajf,Ajf 



(8) 



the conformal dimension of the corresponding CFT a, l 
primary fields. Within the PDT the important functional 
O is well defined for all energy scales and corresponds to 
a much more general paradigm Thus, this connection 
only emerges in the low-energy limit considered here. 

The evaluation of the asymptotic expansion of the 
function on the right-hand side of Eq. Q in- 

volves the use of the general PDT expressions for the cor- 
responding spectral function i3j^(fc, w) of Eq. |(2Jl. After 
some algebra, we find the following general asymptotic 
expansion for the correlation function (O, 



= eilt) J2 
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where Va is the a pseudofermion group velocity at q^^ 
and each term of the X]{a/^} summation is gener- 
ated by transitions from the ground state to the ex- 
cited states with the same values for the current de- 
viations {AJ^,Ajf"}. Note that the asymptotic ex- 
pansion @ has the expected general form, which coin- 
cides with that provided by CFT and used in the stud- 
ies of Refs. El El E m m E3- However, here we 
could obtain an explicit expression for the pre-factors 
Xo — Xoi^Nc, ANs, AJ^ , Ajf) of such an asymptotic 
expansion. It is given by. 



Xo 
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(10) 



where 



Here ^i^' = ^a,a' + Y.,=±li^^)'^aa'iqFa^^lFa') with 

j — 0, 1 involves two two-pseudofermion phase shifts 
whose value is such that in the present low-energy limit 
the square of the shifts in the value of the a, i pseud- 
ofermion Fermi-points given in Eqs. iQ and ((SJ equals 



2tt j 



and the a pseudofermion weight reads. 



(11) 



Aim ^ A^m^^N.^ANs,Aj[,Aj[) = 



)^(gJ-)/2-Q^(g,)/2 + 7^0■-^) 
N* 



n 
n n 



"' n 



,. -2 7r(j - ^) + Qa{qj)/2 
N* 



(12) 



Here N* = Na, N* = + AN^, Qa{qj)/2 is the overaU 
phase shift 10}, and qj e T corresponds to the set of dis- 
crete bare- momentum values in the range qpa.-i < Qj < 



qFa,+i where qFa,L = i-QFa^^lFa.i is the J-ground state 
bare-momentum Fermi point whose deviation is given in 
Eq. (O. The Na dependence of Aa'^^ is such that the 
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quantity given in Eq. l(TT)) is independent of Na- 

Our study refers to the universal part of the asymptotic 
expansion of correlation functions, Eq. lO, and does not 
involve the logarithmic corrections, which are specific to 
each correlation function [l^. The universal character 
of the asymptotic expansion ||^ is such that the value 
of the conformal dimensions only depends on the specific 
correlation function through the values of the four devia- 
tions ATVc, ANg, AJ^, and Ajf of each allowed excita- 
tion J-subspacc. Otherwise, the U, n, and m dependence 
of the two-pseudofermion phase-shift parameters ^, of 
Eq. is specific to the model but is the same for all 
its correlation functions. Importantly, Eqs. H1U|) - (|12|) re- 
veal that the same occurs with the pre-factors xo of the 
asymptotic expansion Indeed, the expression of the 
associated weight of Eq. H12|l involves the overall phase 
shift Qa{q)/2 given in Eqs. Q and 0, which for each 
value of q also depends on the specific correlation func- 
tion through the values of the four deviation numbers 
ANc, ANg, AJf, and Ajf of each allowed J-subspace. 
Moreover, the products of the Aa'^^ expression run in 
the range qpa,-i < qj < qFa,+i, whose limiting values 
deviations are given in Eq. ((SJ and are solely determined 
by the deviations ANa and AJ^ . Thus, the pre-factors 
Xo = XoiANc, ANg, AJ^ , Ajf) value also depends on 
the specific correlation function through the values of the 
four deviation numbers ANc, ANg, AJ^ , and AJf only, 
as the conformal dimensions. Otherwise, the U, n, and 
m dependence of the two-pseudofermion phase shifts in- 
volved in the xo expression is specific to the model but 
is again the same for all its correlation functions. 

In this paper we obtained the pre-factor xo of each 
term of the asymptotic expansion for the correla- 
tion functions of the ID Hubbard model. The form of 
these pre-factors, Eq. pO() . is universal for all correlation 
functions. Their value is controlled by the overall pseud- 
ofermion and hole phase shifts, Eq. Q), through the de- 
pendence on these shifts of the two weig hts aL"'"' of Eq. 
(O and four functionals 2A^ of Eqs. and Con- 
cerning the relation of the latter quantities to the scatter- 
ing mechanisms, note that in the Aq^'"' expression (|12|) 
the bare-momentum products run over the overall phase 
shifts of the a pseudofermion scatterers with bare mo- 
mentum inside the J-ground-state Fermi sea, whose scat- 
tering centers are the c and s pseudofermion and holes 
created at the J-ground-state Fermi points. Furthermore, 
the four conformal dimensions of the GET primary fields 
equal the square of the shifts in the two c and two s 
pseudofermion canonical-momentum Fermi points. The 
four functionals 2AJj and the two weights Aq^'"' also play 
an important role in the finite-energy scattering proper- 
ties, by controlling the unusual spectral properties of the 
model [l,!!^ and real materials 0,01 at all energy scales. 
Thus, our results reveal the connection of the low-energy 
quantities to the scattering mechanisms that control the 



spectral properties at all energy scales. 

While in this paper we considered the ID Hubbard 
model, our results are of general nature for many inte- 
grable interacting problems and therefore have wide ap- 
plicability. Such results provide a broader understan ding 
of existing quasi- ID systems and materials \A lA 14 l2ll|. 
(Recently, the model was used in preliminary investiga- 
tions of the density profiles of ID ultracold fermionic 
atoms in an optical lattice H^.) Indeed, our results re- 
late their low-energy properties to the the general scat- 
tering processes of the objects whose occupancy config- 
urations describe the exotic quantum phases of matter 
corresponding to their different energy scales. This is 
confirmed for finite energies in Refs. [6|,|9(, where the gen- 
eral PDT weight distributions are shown to describe the 
photoemission features of quasi-lD compounds for the 
whole finite-energy band width, whereas the TLL uni- 
versal behavior was observed in quasi-lD materials and 
systems whose low-energy phase is metallic H, ^ 1^ - 
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